Precise electromagnetic tests of ab-initio calculations of light nuclei: States in 10 Be 
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In order to test ab-initio calculations of light nuclei, we have remeasured lifetimes in 10 Be using 
the Doppler Shift Attenuation Method (DSAM) following the 7 Li( 7 Li,a) 10 Be reaction at 8 and 10 
MeV. The new experiments significantly reduce systematic uncertainties in the DSAM technique. 
The r = 2+ state at 3.37 MeV has r = 205±(5) stat ±(7) s „ s fs corresponding to a B(E2[) of 9.2(3) 
e 2 fm 4 in broad agreement with many calculations. The = 2^ state at 5.96 MeV was found to 
have a B(E2[) of 0.11(2) e 2 fm 4 and provides a more discriminating test of nuclear models. New 
Green's Function Monte Carlo (GFMC) calculations for these states and transitions with a number 
of Hamiltonians are also reported and compared to experiment. 



Several ab-initio approaches, i.e., theoretical methods 
that start from "bare" nucleon-nucleon (AW) potentials 
(that reproduce elastic AA scattering data) and em- 
pirical three- nucleon (3 A) potentials, are being devel- 
oped for nuclear physics. Green's Function Monte Carlo 
(GFMC) P,0 and the no-core shell model (NCSM) |,|3 
are two methods that are used most to study p-shell nu- 
clei. They have been successful in reproducing many fea- 
tures, including absolute binding energies and excitation 
spectra^| , charge radii @, 0] , and electromagnetic mo- 
ments [8]. These methods are being developed as suc- 
cessors to the original nuclear shell model [1, [l(| and 
have provided new insight into the origin of the spin- 
orbit force, decoupling of poorly bound neutrons, and 
clustering of nucleons in light nuclei. Several aspects of 
the 3 A force are not yet well constrained, particularly the 
isospin dependence. This is important for the equation 
of state of neutron matter and the properties of neutron 
stars. Thus, testing these approaches by studying some 
of the lightest nuclei can provide insights into some of 
the largest extended nuclear objects in the cosmos. 

The A = 10 nuclei 10 Be, 10 B, and 10 C have long pro- 
vided a stringent test of nuclear models. Kurath 
discussed the need for an unusually strong spin-orbit in- 
teraction to correctly reproduce the sequence of levels in 
10 B. In these nuclei, the p-shell is half-filled and more 
than one state can be created with the same spin-parity 
J"", isospin T, and spatial symmetry. For example, Be 
has two J 7r =2+, T=l states at 3.37 and 5.96 MeV with 
dominant wave-function components 1 D2[442]. These 
states have very different properties and provide a deli- 
cate test of both the Hamiltonians and the calculational 
methods. The Argonne t>i 8 (AV18) AA (TJ interac- 
tion predicts a separation of only 0.23(15) MeV for these 
states, one with a negative [oblate) intrinsic quadrupole 
moment and a small B(E2) decay to the ground state, 
and one with a positive (prolate) quadrupole moment 
which is much more collective. The oblate state is lower 
in energy with just AV18. Calculations including 3A po- 
tentials tend to separate the states and bring the more 
collective prolate-deformed level lower in energy. The 
ratio of electromagnetic decays from these states to the 



ground state is extremely sensitive to structure and varies 
from 4:1 to >50:1 in different nuclear models. In this let- 
ter, we report on new experiments to precisely measure 
lifetimes of bound states in 10 Be. New GFMC calcula- 
tions with a number of realistic forces which explore the 
sensitivity to 3A potentials are also presented. 

To provide significant constraints on the new calcula- 
tions implies achieving <5% accuracy in the E2 matrix 
elements, which is considerably higher than normally re- 
quired for testing nuclear models. The states of interest 
have short lifetimes, a few to hundreds of fs, so Coulomb 
excitation and the Doppler Shift Attenuation Method 
(DSAM) are the appropriate general experimental tech- 
niques. The study of 10 Be following /3-delayed neutron 
decay of 11 Li has recently been shown to also be very 
powerful for this specific case [H, HH • We have concen- 
trated on inferring lifetimes using DSAM, while focusing 
on the reduction of systematic uncertainties associated 
with the method. This technique relies on establishing 
the initial velocity of nuclei in the state of interest and 
then the mean velocity at decay after decelerating in a 
slowing material. From this velocity difference, the life- 
time of the level can be determined if the slowing history 
of the ion is known. Careful selection of the kinematic 
conditions for producing the states of interest at high ve- 
locities where stopping powers are well categorized, con- 
trol of feeding from higher levels, and advancements in 
7-ray detection all are important for moving beyond the 
original >25% measurements from the 1960s [2[l5| and 
improving both precision and accuracy. 

Excited states in 10 Be were populated in the 7 Li( 7 Li, 
a) 10 Bc reaction. Beams of 7 Li ions of up to 50 pnA 
and energies of 8 and 10 MeV were produced by the 
ATLAS accelerator and impinged on thin 7 Li metal and 
7 LiF targets. Target backings of copper and gold of suf- 
ficient thickness to slow the recoiling nuclei to about 
half their original energy were used. The nuclei re- 
coiling along the beam direction were then detected in 
the Argonne Fragment Mass Analyzer (FMA) posi- 
tioned 90 cm downstream of the target and subtending 
1° x 2°. The FMA rejected most non-interacting beam 
particles (> 10 8 suppression) and transported 10 Be ions 
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FIG. 1: (Color online) (Top and middle) Spectra gated on 
10 Be recoil groups with the 7-ray transitions labeled by the 
decaying state. (Bottom) Energy loss versus total energy from 
the ion chamber behind the focal plane of the FMA. The cir- 
cled regions correspond to direct population of distinct 10 Be 
levels. The inset of the top panel shows a level scheme for the 
five bound excited states in 10 Be with 7-ray energies in MeV. 



with charge state g=3 + to the focal plane. The transmit- 
ted ions were identified 50 cm behind the focal plane in a 
20 cm deep, two-electrode ion chamber operated at 15-30 
torr. Gamma rays from the reaction were detected with 
the Gammasphere array 17], consisting of 100 Compton- 
suppressed HPGe detectors in 16 azimuthally-symmetric 
rings mounted at 0=32° to 163° to the beam direction. 
The efficiency of the array (10% for 1-MeV 7 rays) and 
its granularity allowed many systematic tests to be made. 

The positive Q-value (Q=+11.4 MeV) gave 10 Be nu- 
clei high recoil velocities, and thus large Doppler shifts. 
The velocity regime of (3 = v/c = 5-6% is ideal for DSAM 
measurements, as the stopping is 99.99% electronic and 
approaching the "minimum-ionizing" regime for nuclei 
stopping in materials. The current technique only sam- 
pled a small portion of the Bragg curves in the region 
where stopping is best known. The characterization of 
lithium and beryllium ions slowing in lithium, lithium 
fluoride, copper, and gold were taken from the SRIM 



package [18[ , which parameterizes the vast body of stop- 
ping data now available. Modern stopping powers are 
10% higher in this velocity regime than those used in the 
1960s [19]. The initial high recoil energy of ^15 MeV 
is in contrast to the older experiments where it was <1 
MeV and large-angle scattering influenced the direction 
of the emitting nucleus. 

The use of the FMA to select 10 Be nuclei serves sev- 
eral purposes. The most obvious is to provide clean, 
background-free spectra. In addition, due to the two- 
body kinematics of the reaction, a single recoil velocity 
vector is selected where only the 10 Be nuclei traveling at 
0° to the beam axis are studied. Finally, direct popula- 
tion of states can be selected, thereby eliminating cascade 
feeding. In a plot of AE vs E tot for the ion chamber, 
given in Fig. 1, three distinct recoil groups are observed 
corresponding to direct population of the ground state 
(group 1), the 2f state (group 2), and the set of four 
states around 6 MeV (group 3). The wide strip of counts 
below the 10 Be recoils corresponds to 7 Li scattered beam. 
Gating on the direct population of the 2^ state (group 
2), an exceptionally clean spectrum with a single 7 ray 
is observed (Fig. 1, top panel). Selection of the lower 
energy recoils (group 3) gives the decay of the 6 MeV 
states, as presented in the middle panel of Fig. 1. 

For each state there is an optimum target layering, de- 
pending on the lifetime. For maximum sensitivity, equal 
numbers of decays should occur in the target produc- 
tion layer, in the slowing medium, and in the post-target 
flight region. Care was taken to monitor the thicknesses 
of targets and backings, using both weighing and a- 
gauging, before and after the experiments, as this is a 
key source of systematic uncertainty. Several target ma- 
terials and backings were used to investigate systematic 
effects. Lithium metal targets with thickness of 100 - 200 
/ig/cm 2 were produced in an evaporator and transferred 
to the Gammasphere chamber in a vacuum interlock (20j 
maintained at <10 -5 torr. These targets provided the 
simplest 7-ray spectra and allowed the use of high beam 
intensities, as there were few 7 rays produced in contam- 
inant reactions. LiF targets were also used, to provide a 
contrasting stopping power in the target layer. However, 
reactions on the fluorine limited count rates, so these 
data sets were statistically inferior. 

The inference of lifetimes in the DSAM is obtained 
from determining the difference between production and 
emission velocities. If the deacceleration process is 
known, the mean lifetime of the state can then be de- 
duced. The centroid of the 3.37-MeV, 2^ -> 0^ transi- 
tion was measured in each of the 16 angle groups in Gam- 
masphere. The measured centroid as a function of cos(#) 
for a copper target backing is given in Fig. 2(a), com- 
pared to a fit using the relativistic Doppler shift formula. 
Centroids were measured with sufficient precision that 
higher-order terms quadratic in cos((9) were clear, and 
small mechanical mis- alignments could be corrected for. 
The quality of the fit is more clearly illustrated in Fig. 
2(b). Here, the measured centroids are divided by the 
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FIG. 2: (Color online) (a) Measured centroid of the 3.37-MeV 
transition for each of the Gammasphere angle groups com- 
pared to a fit with the relativistic Doppler formula. Lines cor- 
responding to maximum velocity of the 10 Be recoils (dashed) 
and zero velocity (dotted) are shown, (b) Similar to (a) but 
normalized to the best fit function with /3 = 0.05109(19). 



TABLE I: Mean lifetimes from different target and backing 
combinations determined for the 3.37-MeV level in 10 Be. 
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best fit function, a/I - (3 2 / (l-/?cos(0)) with /3=0.05109. 
Although the mean velocity of nuclei at the time of 7- 
ray emission produced in each target-backing combina- 
tion varied widely, they all led to consistent lifetimes. 
The mean lifetimes extracted from the analysis of the 
different backings are summarized in Table I. The cal- 
culation of the velocity-lifetime relationship was made 
by sub-dividing both the production target and backing 
into 30 layers and determining the average velocity pro- 
file and the probability of emission in each layer. The 
inferred lifetime did not change when a finer layering 
was used. Both the statistical uncertainties (from the 
fit to the Doppler shifts) and the systematic uncertain- 
ties (from thickness and stopping power effects) varied 
considerably for the different target and backing arrange- 
ments. However, the many varied measurements allowed 
appraisal of systematic uncertainties. 

The 3.37-MeV level was measured in seven separate ex- 



periments. The weighted mean value is t— 205±(5) s t a t 
±(7) S y S t fs which gives a reduced transition probability 
B(E2i)=9.2(3) e 2 fm 4 . Our lifetime is about 10% longer 
than the mean of previous experiments plj which span 
the range 110 fs to 260 fs at the la level. The leading con- 
tribution to systematic uncertainty was the target layer 
thickness, especially for the lithium metal targets. In the 
future, Rutherford backscattering of the beam will be 
monitored during the experiment to determine the tar- 
get thickness in live time at the point where the beam 
hits. The next leading term was the uncertainty in the 
stopping power of the slowing medium, which is now at 
the < 2.5% level for this velocity regime [l8l ]. 

Four bound states near 6 MeV were all populated and 
spectroscopic information was extracted. For this discus- 
sion, we focus on the second J 7r =2 + state at 5.96 MeV. 
This level is part of a close- lying doublet with a J^l - 
level that is only ^1.5 keV higher and the 7 decays from 
the two states could not be resolved. Thus, analysis 
of the 2.59-MeV and 5.96-MeV transitions yielded dif- 
ferent "effective" lifetimes of 50(8) fs and 30(8) fs, re- 
spectively, highlighting the doublet nature of the transi- 
tions. Unfolding the doublets required knowledge of the 
7-branching ratios from each level and their relative pop- 
ulation intensities. Combining these two quantities and 
the measured "effective" lifetimes, the second 2 + state 
was found to have r= 59±(15) fs. For the 2\ — ► Of 
branch, we used the precise value of 6.1(1)% measured in 
neutron capture (22j. In the decay to the ground state, 
this gives £(£2j)=0.11(2) e 2 fm 4 . Complete details of 
the analysis will be given in a future paper [23| . 

In parallel to the experimental work, extensive GFMC 
calculations of states and transitions in 10 Be were carried 
out using AV18 and a variety of 3iV forces, see Table II. 
The absolute GFMC energies are believed to be accurate 
to 1-2% for a given Hamiltonian, with somewhat larger 
errors for the excitation differences. The GFMC method 
for evaluating electromagnetic decays between states is 
detailed in Ref. [8j. Calculations with the AV18 inter- 
action alone predict near-degeneracy for the two J 7r =2 + 
states, with the oblate state slightly lower in energy. To 
explore the sensitivity to 3N potentials, calculations were 
performed with AV18 plus Urbana-IX (UIX), Illinois-2 
(IL2), and Illinois- 7 (IL7) Inclusion of the 3iV forces 
lifts the degeneracy between the 2 + states, slightly for 
UIX, more so for IL2 and approaching the experimen- 
tally observed spacing for the IL7 potential. 

Similar trends were found for the B(E2) transition 
strengths to the ground state. With just the AV18 
potential, the calculations overpredict both transition 
strengths. UIX provides a good description of the de- 
cay of the 2% state, yet significantly overpredicts the de- 
cay of the 2f state. With the IL potentials, it appears 
that as the 2 + states become farther apart in energy, 
the predicted B{E2[) strengths tend to approach the ex- 
perimental values. Considering both the energies and 
transition strengths, the IL7 potential provides the best 
description. However, while the summed E2 transition 
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TABLE II: GFMC calculations of the 10 Be ground-state energy E gs , excitation energies E x (both in MeV), and B{E2[) tran- 
sitions in e 2 fm 4 .The calculations were done with a variety of potentials to explore the sensitivity in predicting electromagnetic 
matrix elements. 



H AV18 AV18+UIX AV18+IL2 AV18+IL7 Expt. 

\E ga (Q + )\ 50.1(2) 59.5(3) 66.4(4) 64.3(2) 6T98 

E x {2+) 2.9(2) 3.5(3) 5.0(4) 3.8(2) 3.37 

E x (2+) 2.7(2) 3.8(3) 5.8(4) 5.5(2) 5.96 

B(£2;2+ ^ 0+) 10.5(3) 17.9(5) 8.1(3) 8.8(2) 9.2(3) 

B(£2;2+^0+) 3.3(2) 0.35(5) 3.3(2) 1.7(1) 0.11(2) 

Y,B{E2) 13.8(4) 18.2(6) 11.4(4) 10.5(3) 9.3(3) 



strength from the 2 + states to the ground state predicted 
by the IL7 potential is in very good agreement with ex- 
periment, the distribution of strength is not correct, with 
too much strength computed for the upper state. This 
shortcoming is still being investigated. It could possibly 
be remedied by further refinement of the Illinois three- 
body potentials, or by including more extended 2 + com- 
ponents arising from particles in the sd-shcll. 

A wide variety of other models have also been used 
to investigate 10 Be. In the traditional p-shell model 
with unmixed A"-quantum number, Millener [25[ predicts 
B(E2;2+ -> 0^)=11.8 e 2 fm 4 and the ratio of B(E2)s as 
50:1. Both the NCSM with the CD-BONN, AW poten- 
tial and the Microscopic Cluster Model (MCM) [2f3] 
tend to underpredict the B(E2; 2^ — > Of) decay, with 
NCSM predicting 6.6 e 2 fm 4 and the MCM predicting 6.1 
e 2 fm 4 . However, both predict a weak B(E2; 2\ — » fT*~) 
decay, ~0.13 e 2 fm 4 , in good agreement with experiment. 
Thus, at present, no model provides a perfect descrip- 
tion of 10 Be and reconciling the predictions from these 
different models is an ongoing challenge. 



In conclusion, the Doppler Shift Attenuation Method 
(DSAM) has been refined to precisely determine the life- 
time of the first excited J T =2 + state in 10 Be. Care was 
taken to control the systematic uncertainties usually as- 
sociated with DSAM. High precision lifetime measure- 
ments now appear feasible in not only light nuclei, but 
also a range of systems which can be produced through 
two-body kinematic reactions. The lifetime of the 2^ 
state in 10 Be was also measured. New ab-initio calcu- 
lations show the sensitivity of the transition matrix el- 
ements to nuclear structure, especially to the form of 
three-body forces. In GFMC calculations the overall 
B(E2) decay strength to the ground state can be well 
reproduced in calculations, but its distribution is still 
difficult to account for accurately. 
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tex computer of the Argonne Mathematics and Computer 
Science Division. 
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